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For the initial fields of the density contrast and peculiar velocity, we theoretically calcu- 
late the 'differential' and 'integral' length scales, i.e., statistical measures that respectively 
characterize the small- and large-scale fluctuations of a random field. These length scales 
and the associated mass scales explain the length and mass scales observed for (1) halos of 
young galaxies at z > 5, (2) halos of galaxies at 2 ~ 0, and (3) the largest structures in the 
galaxy distribution at 2: ~ 0. We thereby discuss that such observed scales are fossil imprints 
of the characteristic scales of the initial fields. 

Subject Index: 464, 469 

§1. Introduction 

Since structures of the Universe are formed from gravitational instability of the 
initial fluctuations, some characteristics of the structures are traced back to those 
of the initial fluctuations. This should be especially the case for the characteristics 
determined in the linear or quasi-linear stage of the structure formation. Both the 
stages are described with the initial peculiar velocity v by the first-order Lagrangian 
perturbation.lll'''^ The linear stage is also described with the initial density contrast 
5 by the first-order Eulerian perturbation.'^''3J 

The characteristics of the initial fields of the density contrast 5 and peculiar 
velocity v are studied here. Since these fields are random, we use statistical measures 
that have been developed to characterize turbulent random fields.® Specifically, for 
the first time, the 'differential' and 'integral' length scales are used to characterize 
the initial fields. We calculate the length scales and associated mass scales. They are 
compared with the scales observed for galaxies and their large-scale distribution, in 
order to explore the fossil imprints of the initial fluctuations among these structures. 

§2. Basic framework 

2.1. Cosmological parameters 

For the cosmological parameters, we assume their standard values: the Hubble 
constant Hq = 70kms~^Mpc~^ {h = 0.7), the matter density Q^a = 0.14/i~^, the 
baryon density i7b = 0.023/i~^, and the temperature of the cosmic microwave back- 
ground TcMB = 2.73K.ISJ' The matter density Q^a is dominated by the cold dark 
matter so that the structure formation proceeds from small to large scales. 
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2.2. Density and velocity fluctuations 

The density contrast 5{x,t) at comoving position x and at time t is defined as 
6{x,t) = p{x,t)/{p) — 1 with the mass density p{x,t). Here, (•) denotes an average 
and {p) = ZQ^Hq /SttG. During the hnear stage, the Unear growth factor D{t) 
describes the evolution as 5{x,t) = D{t)6{x). 

The density contrast d{x) at the initial time tin makes up a random field, which 
is related to its Fourier transform d{k) as 

/dh 
5ik)eMik-^)j^- (2-1) 

Since 5{x) is homogeneous, isotropic, and Gaussian, at least as a good approxi- 
mation, its statistics are determined by the initial power spectrum Ps{k) for wave 
number k = \k\: 

{6{k)6{k'y) = {27:f6B{k - k')P5{k). (2-2) 

Here, S{k)* is the complex conjugate of S{k), and 6D{k) is Dirac's delta function. 

For the growing mode of the fluctuations, the peculiar velocity v{x) at t = tin 
is defined by using 6{k): 

v{x) = D{t,^) I exp(zfc . (2-3) 

(Peebles,'^ p. 514), where D{t) is the time derivative oi D{t). The factor k/k'^ implies 
that v{x) is dominated by larger scales than 6{x). Length scales for v{x) are larger 
than those for S{x) (see 

The peculiar gravitational potential cj){x) at t = tin is also defined by using 6{k): 

(Peebles,'^ p. 114), where z{t) is the redshift. We have v{x) oc — V0(a;). That is, 
at t = tin, the peculiar velocity v(x) for the growing mode is proportional to the 
peculiar gravitational acceleration —'V(j){x). 

Usually, the initial fields are smoothed with a window function.'^''^' The reason 
is the significant small-scale fluctuations retained by the cold dark matter. For 
example, Mouri and TaniguchP* smoothed 6{x) and v{x) to calculate some of their 
characteristic scales. However, any smoothing imposes its scale as an additional 
artificial scale. The present study is to be among the first to calculate scales of the 
unsmoothed initial fields. 

2.3. Perturbation theories 

The first-order Eulerian perturbation 6{x,t) = D(t)6{x) has been mentioned in 
^2.21 This holds for regions in the linear stage. The absolute value of S{x, t) averaged 
over such a region is <C 1. 

The first-order Lagrangian perturbation, i.e., the Zel'dovich approximation,!^'' 
describes the position x(t) of a matter element with its initial position xi^ and its 
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initial peculiar velocity v^xi^): 

D(t) , , , , 

x{t) = Xin + j-^v{xi^). (2-5) 

This holds for regions in the linear or quasi- linear stage. The absolute value of 6{x, t) 
averaged over such a region is < 1. In other words, Eq. (j2-5p is valid so far as the 
matter elements from different initial positions do not come across one anotherP'^^ 

2.4. Models for structure formation 

The perturbation theories allow us to relate scales of the initial fields to scales of 
structures determined in the linear or quasi-linear stage of their formation. However, 
these structures eventually become nonlinear. To such self-gravitating objects with 
6{x,t) » 1, the perturbed initial fields are extended by using models. 

For scales of the density field, we use the peak model of Kaiser.l2I''[3i'Eli This is 
based on the Eulerian perturbation 5{x,t) = D{t)6{x). When D(t)6{x) averaged 
over a region exceeds some critical value of order unity, the region is assumed to 
collapse and form a self-gravitating object. It follows that the structure formation 
occurs preferentially around peaks of the initial density field. 

For scales of the velocity field, we use the adhesion model of Gurbatov et ah'^'-EJ 
This is based on the Lagrangian perturbation in Eq. (|2-5p . When matter elements 
from different initial positions come across one another, the adhesion model assumes 
that the matter elements adhere together, mimicking the gravitational interaction 
and the resultant formation of a self- gravitating object. 

Thus, for scales of the density and velocity fields, we use different models, al- 
though they represent the same dynamics. Details of the models are not important 
here because we only require their qualitative predictions. 



§3. Characteristic scales 



3.1. Definition and meaning 

Consider a homogeneous isotropic random field f{x) with {f{x)) = 0. Its two- 
point correlation at scale r = \r\ is {f{x)f{x + r)), which is used to define the 
differential and integral length scales. Since these length scales possess not only the 
dimensions of length but also particular meanings, they deserve to be compared with 
the observed scales. 

The differential length scale A is defined by using the correlation and its second- 
order derivative at r = as 

1/2 



A = 7r 



d^.{f{x)f{x + r))\r=o 



(3-1) 



This is identical to A = ^ [{f{xf)/{{dxff)Y^'^, where X is any one-dimensional cut 
of the space x (Monin and Yaglom,® pp. 25 and 35). The original definition lacks 
the factor vr. Our definition with the factor vr is essential so that A has a particular 
meaning. That is, if f{x) is Gaussian as in the case of the initial density and velocity 
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fields, A is exactly identical to the mean interval between successive zero crossings 
of f{x) along any one-dimensional cut of the spaceP^ Although past studies used 
scales that differ from A only by numerical factors P^'EJ-CS difference is essential 
so that A is a new statistical measure with the new particular meaning. 

The integral length scale L is defined by integrating the correlation over all r as 

^ J,°"{f{x)f{x + r))dr 
" ifW) ^ 

(Monin and Yaglomj^ p. 34). This is the characteristic length scale of significant 
correlation. While {f{x)f{x-\-r))/{f{x)'^) is significant at r < L, it is not significant 
at r > L. An elementary example is a correlation that is proportional to an expo- 
nential function exp(— r/L) or a step function 9{L — r). The situation is the same for 
correlations in the initial density and velocity fields (see To characterize such 
correlations, past studies used scales oc A.'^''^ However, they are not the length 
scales of significant correlation because A is defined only at r = (Eq. ()3-ip ). In 
fact, we usually have X <^ L. 

3.2. Condition for existence 

There is a condition required for the existence of the length scales A and L. They 
have been defined with the correlation function, which is in turn defined with the 
power spectrum. For the initial density and velocity fields studied here, the length 
scales are in the forms of 



A oc 



j,°^k-^Ps{k)dk f j^k^'^m)dk 

/o°°A:-+2p,(A;)dfcJ ^""'^ J^kmPs{k)dk ' ^ ^ ^ 



with integer m (Eqs. (|4-2p . (j4-4p . (j4-6p . and (|4-7p ). The above integrals are required 
to converge. This condition is important especially because we use the unsmoothed 
initial fields f ^O) . 

To describe the initial power spectrum, we have Psik) cx T{k)'^k'^''. Since the 
transfer function T{k) is dominated by the cold dark matter, it has the asymptote^^ 

limr(/t) = l and lim T(/t) oc i^^. (3-4) 

Here, we have ignored the cutoff at a very high k that reflects the particle physics 
of the cold dark matterjl^Sli'lIII for which little is known. The hot dark matter and 
baryon are not dominant and hence affect the asymptotes only through the prefactor 
of ln{k)/k'^^^'^^ On the other hand, for the primordial power spectrum k^" , the 
standard model is the Harrison-Zel'dovich spectrum with Ug = 1. However, the 
actual spectrum for the observable Universe is tiltedP 

ns = l- As with As ~ 0.04, (3-5) 



because of some deviation from exact exponential expansion of the Universe during 
the corresponding stage of the inflation. 
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Table I. Length scales L and A and associated mass scales 
M{L) = A'k{p)L^/?, and M(A) = 47r(p)AV3 for = 
0.96 ± 0.01, nmh"^ = 0.14, n^h'^ = 0.023, and Tcmb = 
2.73 K. The scales are smaller for the greater value of Ug. 



Length Mass Eq. 



(Mpc) [Mq) 



Ls 




M{Ls) 


2 0-° * 


X 10" 


(|T2)l 


A„|| 








X 10" 


(l4-6a|) 


Au^ 


2 S-O-^ 

^•■^+0.3 




2 fi-"'^ 

^•"+0.9 


X 10^2 


(IX^ 


Ai; 


2 0-°-2 

^•"+0.1 


M{Xv) 


1 2~"'^ 

-^•^+0.4 


X 10^^ 


(|4-6b|) 


Lvj^ 


130.;^; 




o (5-0.2 


X 10^^ 


(IX^ 


Lv 


87.;?: 


M{Lv) 


^■^+0.0 


X 10" 


(|4-7b|l 



Hence, in Eq. p-3p . while the integrals for A converge if m is —1 or 0, those for 
L converge if m is 0, 1, or 2. This condition is satisfied by the density integral scale 
in Eq. ()4-2p . velocity differential scale in Eq. ()4-6p . and velocity integral scale in Eq. 
(|4-7p , but not by the density differential scale in Eq. ()4-4p . The last scale does not 
existol It should also be noted that the condition is not satisfied by either A or L 
for the initial field of gravitational potential. 

§4. Results and discussion 

Table |I] lists the numerical values of the length scales L and A as well as of the 
associated mass scales M{L) = 47r(p)L^/3 and M(A) = 47r(p)A^/3. 

For the primordial power spectrum k"'" , since the Ug value is crucial to the 
convergence of the integrals in Eq. (|3-3p , we have used the value with an uncertainty, 
Hg = 0.96 lb COl-l^ Table U shows that the corresponding uncertainties are relatively 
less significant in the larger scales. 

For the transfer function T{k), we have used the fitting formulae of Eisenstein 
and Hu.'^ Their input parameters are f2^h'^, Q^^h?^ and TcmB; for which we have 
used the standard values in ^2.lF*)l Since Q\yh? <^ ^mh?, features of the baryon 
acoustic oscillation^I^ are not so strong in T{k). Then, the overall shape of Ps{k) oc 
T{k)'^k^'' determines the values of the scales. The hot dark matter, i.e., neutrino, 
has been safely ignored because its particle mass, albeit not exactly known, is small 
enoug 

Hereafter, the individual scales are studied along the structure formation from 
small to large scales, or equivalently, from early to late times. 



*' If we assume the Harrison-Zel'dovich model spectrum with Us — 1, only the velocity integral 
scale exists. For the other scales, the integrals do not converge because small-scale fluctuations of 
the unsmoothed initial fields are too significant (see also t|2.2|l . Thus, it is essential to use the actual 
spectrum with = 1 — As — 0.96. 

**•* If we assume an uncertainty of ±0.01 in Omh^ , ±0.001 in Obh^ , or ±0.01 K in Tcmb, it leads 
to an uncertainty of about ^20%, ±3%, or ±2% in each of the mass scales. 
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Fig. 1. Two-point correlations {S{x)S{x + r)) , {v{x) ■v{x + r)} , {v^^{x)v^^{x + r)) , and {v±{x)v±{x + 
r)) at the initial time tin as a function of scale r. They are normalized by the values at r = 0. 
The arrows indicate the integral length scales Ls and Lv- 



4.1. Density integral scale 

By using Eqs. (|2-ip and ()2-2p . from the initial power spectrum Ps{k), we obtain 
the two-point correlation of the initial density contrast 6{x) as 



{5{x)5{x + r)) = 4tt k'^jo{kr)Ps{k)dk (4-1) 
Jo 

(Peeblespl' p. 509), where ^0(2:^) = sin(x)/x is a spherical Bessel function. The profile 
of {5{x)6{x + r)) / {5{x)'^) is shown in Fig. [TJ 

Equations (|3-2p and (|4-ip are used to obtain the integral length scale for the 
initial density contrast: 

_ ^j^kpmdk 

^ 2j;^k^Ps{k)dk ^ ' 

Table |l] shows that the length scale is about 0.02 Mpc and the associated mass 
scale M{Ls) is about 2 x 10^ Mq. We regard them as lower limits. The actual values 
could be much greater. This is because we have ignored the high-fc cutoff in Ps{k) 
( ^3.2p . To the other scales in Table [H the cutoff is not important!*^ 

The mass scale M(Ls) serves as the typical mass for overdense regions around 
high density peaks of the initial fluctuations. If such a peak at Xp has the density 
contrast S{xp), the mean density contrast at distance r is 



{6{x)6{x + r)) 
7^ 



{S{xp + r))= '"^,Z':j/^' 6ixp) (4.3) 



*-* If we assume the cutoff exp[— (fc/fec)^] at a relatively low wave number kc = 10'' Mpc" , 
M{Ls) is 37 times greater, while M{Xv) and M{Xv) are 6.1 times greater. If kc = lO'^ or 10* Mpc~^, 
M(Ls) is 14 or 5.8 times greater, while M(Ai,) and M{\v) are 3.7 or 2.4 times greater. For each 
kc, M{Lv) and AI{Lv) remain the same. 
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(Kaiser,-** see also Peebles,^' p. 623). This is identical to the radial profile of the 
mean overdensity. Figure [1] shows that {6{x)6{x + r)) / {6{x)'^) is significant up to 
r Ls The mean overdensity {S{xp + r)) is thereby significant up to the 

radius of r ~ Ls- Then, Ls serves as the typical radius for the overdense regions, 
and hence M{Ls) = i7r{p)L^/3 serves as the typical mass. 

In practice, any overdense region consists of many density peaks. Their typical 
size is just the differential length scale for 6{x), i.e., mean interval between its zeros 
along any one-dimensional cut of the space ( ^3.ip : 

•^5 = ^ TOO ,^ „ X „ ■ (4-4) 



f^k^Ps{k)dk 

Since the integral k^Ps{k)dk diverges, the scale is always zero. The sizes of 
those density peaks are very small. 

The overdense region around a density peak S{xp) collapses and forms a self- 
gravitating object, according to the peak model in ^2.4[ If S{xp) is high enough, 
the collapse occurs early enough so as not to be affected by other overdense regions. 
The object mass is typically M{Ls) that has a value at > 2 x 10^ Mq. Being more 
massive than 'minihalos' with 10^-10^ Mq of the first stars ,'^''^31 this class of objects 
formed at early t is considered to be a young stage of galactic halos once existent at 
high z. 

The representative young galaxies observed at z > 5 are 'Lya emitters' that are 
characterized by intense Lya emission from ionized gas around hot massive stars. For 
the individual Lya emitters observed so far, the range of the stellar mass is down to 
10*^ Mq.IEJ'ES This is consistent with the total halo mass of M (Ls) at > 2 x 10*^ Mq, 
even though is several times greater than f?b (El]) and a fraction of baryon is 
in the gas. 

4.2. Velocity differential scale 

By using Eqs. (|2-2p and ()2-3p . from the initial power spectrum Ps{k), we obtain 
the two-point correlation of the component f y of the initial peculiar velocity v(x) 
that is parallel to the separation r: 

2ji{kr) 



{v\\{x)v\\{x + r)) =4TTD{tinf / 

Jo 



kr 



Ps{k)dk. (4-5a) 



Here, ji (x) = sin(x) /x'^ —cos{x) jx is a spherical Bessel function. Also for the velocity 
vector v{x), the two-point correlation is obtained as 

{v{x) ■ v{x + r)) = 47rl)(tin)2 / jo{kr)Ps{k)dk (4-5b) 

Jo 

(G6rski,l2D see also Monin and Yaglom,^!' p. 51). We show the profiles of (wy {x)v\\ {x+ 
r)) / {v\\{x)'^) and {v{x) ■ v{x + r)) /{\v{x)\'^) in Fig. [TJ 

Equations ()3-ip and (j4-5ap are used to obtain the differential length scale for the 
v\\ component: 



_5£^Psik)dk_-^'^' 
3f^k^Psik)dk 



(4-6a) 
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Table H] shows that the length scale A^|^ is about 2 Mpc and the associated mass scale 
M(A^ll) is about 5 x 10^^ Mq. 

For reference, we use Eqs. (j3-ip and ()4-5bp to obtain the differential length scale 
for the velocity vector v: 

1/2 

(4-6b) 



A 



V 



vr 



k^Ps{k)dk 



The meaning of A explained in ^3.11 does not hold if A is for a vector. Hence, A^; is 
not studied here, albeit mentioned in ^4.31 (see also Appendix |A|) . 

The mass scale M(A^||) serves as the typical mass for clustering toward local 
minima of the initial gravitational potential (/)(a;). This clustering is described with 
the initial peculiar velocity v by the Lagrangian perturbation in Eq. ()2-5p . Figure [2] 
shows the initial fluctuations along a one-dimensional cut of the space. As indicated 
by arrows, there are clustering motions according to the v^^ component that is parallel 
to the one-dimensional cut. Zeros of the v^^ component are centers or boundaries of 
the clustering. They correspond to the local minima or maxima of (pix) because 
v{x) oc — V0(a;). The mean interval between those zeros is A^n ( ^3.ip . which serves 
as the typical clustering radius. Hence, M{Xy^^) = 47r(/j)A^|| /3 serves as the typical 
mass that undergoes the clustering. 

To be exact, the one-dimensional cut does not necessarily pass through the 
clustering centers in the three-dimensional space. This fact is nevertheless not serious 
at all, as demonstrated in Appendix [Bl 

The matter elements that have clustered around a local minimum of the initial 
gravitational potential (j){x) interact with one another. They eventually form a self- 
gravitating object. This is just the prediction of the adhesion model in ^2.41 (Sahni 



density contrast 




gravitational potential 



10 15 
position (Mpc) 
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Fig. 2. Density contrast 5{x), peculiar velocity v\i{x), and peculiar gravitational potential 0(a3) at 
the initial time tin along a one- dimensional cut of the space. They are in arbitrary units. The 
arrows indicate the clustering motions. We also show the difTerential length scale X^,. . 
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and Coles, ^' see also Kofman et al.^^'). Such a self-gravitating object has the typical 
mass M(A^|| ). 

We have studied young galaxies formed at high z ( ^4.1|) . Their halos with mass 
~ M(Ls) represent the matter elements that cluster to form the objects with mass 
~ M(A„ll) 3> M{Ls) (Table [J). On the other hand, it is widely considered that the 
halos of young galaxies accrete matter and merge with one another to form the larger 
halos of galaxies observed at z ~ Op^ They are identical to the objects with mass 
~ M(A^ll) if the accretion and mergers occur within each potential well of (/'(a;). 

The galactic halos observed at z ~ had been formed almost completely at 
z > 0. Most of these halos have been accreted into the larger halos of galaxy groups 
or clusters but have survived as their identifiable substructures P^'EU There also 
remain isolated galaxies that do not belong to groups or clusters. 

The representative galaxies observed at z ~ are those at around the Schechter 
luminosity ~ W^^Cq. Although the less luminous galaxies are more numerous, 
galaxies at around dominate the luminosity integrated over all galaxies (Peebles,'^ 
p. 120), judging from the observed shape of the luminosity function. Since the mass- 
to-light ratio is not sensitive to the galaxy luminosity, it follows that the galaxies 
at around £* share most of the mass. From the observed gravitational lensing and 
motions of satellite galaxies, it has been established that each isolated galaxy at 
has the total halo mass of (1-2) x 10^^ MqW^^^^ The same halo mass is obtained 
for galaxies at in groups and clusters, except for those in the cluster core where 
the tidal field is too strong.l^S'l^l They are close to M{X^^^) ~ 5 x 10^^ Mq. 

4.3. Velocity integral scale 

Equations ()3-2p and ()4-5|) are used to obtain the integral length scales for the 
initial peculiar velocity v{x) from the initial power spectrum Ps{k): 

L^ll = 0, (4-7a) 
_ n k-^ Ps{k)dk 

- 2f-Ps{k)dk ■ ^'-'^^ 

The length scale is always zero because {v\\{x)v\\{x + r)) changes its sign (Fig. 

d]), while the length scale Lv is about 90Mpc (Table |Il)E3 

Table |T] shows ^ At) ~ 2 Mpc. Beneath the small-scale fluctuations that 
dominate the differential length scale Xv (see Fig. [2]), there is hidden the large-scale 
correlation that dominates the integral length scale L-d (see Fig. [T]). 

The length scale serves as the typical size for the largest structures of the 
Universe such as those in the galaxy distribution. Even at z ~ 0, their formation 
is ongoing and in the quasi-linear stage. The Lagrangian perturbation in Eq. (|2-5|) 
implies that the formation is due to large-scale correlated motions in the initial 
velocity field v{x)^^^ As shown in Fig. [H the correlation of v{x) is significant up 
to the scale r ~ L-d ( ^S.ip . The structure formation is accordingly significant. Above 
the scale r ~ Lv, the correlation is not significant, and hence the structure formation 



*-* The mass scale M{Lv) ~ 47r(p)L|,/3 is not considered here because the corresponding 
structures are not self-gravitating even at z ~ (see below). 
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is not significant. Thus, on length scales for the significant structure formation in 
the quasi-linear stage, there is an upper limit Lv- 

We underline that the above discussion is statistical. Some structures could be 
formed far above the scale r ~ L^, where {v{x)-v{x + r))/{\v{x)\'^) is not significant 
but is not yet absent (Fig. [T|). The initial velocity field v{x) could be correlated over 
very large scales in some exceptional regions, which could form very large structures. 
However, such exceptional structures are not statistically important. 

The galaxy distribution observed at 2 ~ has large-scale structures such as 
filaments, walls, and voids. Statistical signatures for the existence of the structures 
are significant up to 50-100 Mpc, above which the galaxy distribution is almost 
uniform .ED"t2Sli This length scale for transition to uniformity is close to Lv — 90 Mpc. 
Although structures much larger than are known, i.e., 'great walls','^ they are 
exceptional and are not significant to the statistics.!^ 

The observed galaxy distribution is also used to reconstruct the peculiar velocity 
field at z ~ 0,E3'I3SJ where we see large-scale correlated motions from voids toward 
filaments and walls. The length scales of these motions, albeit not explicitly given 
in the literature, are up to about Ld ~ 90 Mpc. 

According to the adhesion model in ^2.41 the structure formation has two stages 
(Sahni and Coles,'21' see also Kofman et al.^^). The first stage is clustering of matter 
elements toward local minima of the initial gravitational potential 4'{x) and the 
resultant formation of self-gravitating objects. The second stage is clustering of 
these objects due to a large-scale correlation in the initial gravitational acceleration 
—'V(j){x). Up to some characteristic scale of the correlation, the structure formation 
is significant. By characterizing the initial peculiar velocity v(x) oc —'V(j){x) with 
the new scales and L-y, we have found that, while formed in the first stage are 
halos of galaxies with mass ~ M(At,||) ( ^4.2p . formed in the second stage are the 
largest structures in the galaxy distribution with size ~ L^- 

§5. Conclusion 

For the initial fields of the density contrast 5 and peculiar velocity v, we have 
calculated new statistical measures: the differential length scale A, the integral length 
scale L, and the associated mass scales M(A) and M{L). Through the first-order 
Eulerian or Lagrangian perturbation, these scales are related to scales for the struc- 
ture formation. We have found that they explain the scales observed for galaxies 
and their distribution. 

The mass scale M{Ls) at > 2 x 10® Mq serves as the typical mass for overdense 
regions around high density peaks of the initial fluctuations, which explains the 
typical mass for halos of young galaxies observed at z > 5. 

The mass scale M{Xy^^ ) ~ 5 x 10^"*^ Mq serves as the typical mass for clustering 
toward local minima of the initial gravitational potential, which explains the typical 
mass for halos of galaxies observed at z ~ 0. 

The length scale — 90 Mpc serves as the typical size for the largest structures 
formed in the Universe, which explains the length scale for transition to uniformity 
in the galaxy distribution observed at z ~ 0. 
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Therefore, these observed scales are fossil imprints of the characteristic scales of 
the initial fields. While the halos of young galaxies with mass ~ M{Ls) are extinct 
at z ~ 0, the halos of galaxies with mass ~ M{Xy^^) and the largest structures with 
size ~ Ld are existent and observable at z ~ 0. 

Appendix A 

Other Characteristic Scales of the Peculiar Velocity 



The differential and integral scales also exist for the v± component of the initial 
peculiar velocity that is perpendicular to the separation r. Although this v± compo- 
nent is not directly relevant to the structure formation studied here, we summarize 
its scales for reference. 

By using Eqs. (|2-2p and (|2-3p . we obtain the two-point correlation of the v± 
component^^' 

{v^{x)v^{x + r)) = A7rD{t,^f ^2^p,(k)dk. (A-1) 

The profile of {v±{x)v±{x + r)) / {v±{x)'^) is shown in Fig. [TJ Equations (|c{-l|) . 
and (|A-1|1 are used to obtain the length scales 

5rps{k)dk^'/' 



Xy=TT 



k^Ps{k)dk 
Stt k-^ Ps{k)dk 



(A-2) 



Table [J lists the numerical values of these length scales and of the associated mass 
scales. 

Appendix B 

Clustering Mass Obtained from the Gravitational Potential 



We have discussed the mass scale M{Xy^^ ) as the typical mass for clustering 
toward local minima of the initial gravitational potential (/>(£c), by ignoring the fact 
that a one-dimensional cut does not necessarily pass through the local minima in the 
three-dimensional space ( ^4.2p . Here is a more exact discussion that leads to almost 
the same result. 

We calculate the number density of the clustering centers, i.e., number den- 
sity nrnin ,^ of local minima of the initial gravitational potential (jy^x). For a three- 
dimensional random scalar field that is homogeneous, isotropic, and Gaussian as in 
the case of 4>{x), the number density nmax of the local maxima was calculated by 
Bardeen et al^ To their Eqs. (4-6) and (4-11), we substitute the power spectrum 
of 4>{x) obtained from Eqs. ()2-ip . ()2-2p . and ()2-4p . which is proportional to Ps{k)/k^. 
Then, 



29 - qVg 



3J-Ps{k)dk f' 
J^k^Ps{k)dk\ ■ ^ ^ 



12 



H. Mouri and Y. Taniguchi 



We have oc (Eq. ()4-6ap ). although does not have any particular meaning. 
The number density nmm,!/) leads to the typical mass {p) / fiyain^cj) for clustering toward 
local minima of the initial gravitational potential ipix), which is only 15% greater 
than M(At,||) = 47r(p)A^ /3. In fact, our cosmological parameters yield {p) / n^ia,(f> = 



5.7"};^ X 10^1 Mq, while M{X,, ) = S.OTj;^ x 10^^ Afg 



1.3 
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